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With the advent of neoclassically optimised stellarators, optimising stellarators for turbulent transport is an important 
next step. The reduction of ion-temperature-gradient-driven turbulence has been achieved via shaping of the magnetic 
field, and the reduction of trapped-electron mode (TEM) turbulence is adressed in the present paper. Recent analytical 
and numerical findings suggest TEMs are stabilised when a large fraction of trapped particles experiences favourable 
bounce-averaged curvature. This is the case for example in Wendelstein 7-X [C.D. Beidler et al Eusion Technology 17 , 
148 (1990)] and other Helias-type stellarators. Using this knowledge, a proxy function was designed to estimate the 
TEM dynamics, allowing optimal configurations for TEM stability to be determined with the STELLOPT [D.A. Spong 
et al Nucl. Fusion 41 , 711 (2001)] code without extensive turbulence simulations. A first proof-of-principle optimised 
equilibrium stemming from the TEM-dominated stellarator experiment HSX [F.S.B. Andersonef al. Fusion Technol. 
27 , 273 (1995)] is presented for which a reduction of the linear growth rates is achieved over a broad range of the 
operational parameter space. As an important consequence of this property, the turbulent heat flux levels are reduced 
compared with the initial configuration. 


I. INTRODUCTION 

The development of stellarators has taken great strides since 
their inception by Lyman Spitzer Jr. in 1951—. Stellara¬ 
tors are inherently 3-dimensional, so that the configuration 
space of possible equilibria is very large. Different tech¬ 
niques of optimisation have been employed to And equlibria 
with certain desired features within this configuration space. 
The reduction of neoclassical transport down to levels of 
tokamaks for example has been achieved through the intro¬ 
duction of quasi-symmetries^i^ or variants of omnigeneityi^, 
which made stellarators competitive with tokamaks regard¬ 
ing the expected levels of transport. In addition, other stel¬ 
larator design features are being addressed via optimisation, 
too, such as the confinement of fast particles^. In neoclas¬ 
sically optimised stellarators such as Wendelstein 7-X (W7- 
X)^ or the quasi-symmetric stellarators HSX (Helically Sym¬ 
metric Experiment)^ and NCSX (National Compact Stellara¬ 
tor Experiment)i2, the turbulent transport is expected to be the 
dominant transport channel in a large part of the plasma, just 
as in tokamaks. While NCSX has not been built and W7-X 
is only about to start operation, HSX is already running and 
has shown that the neoclassical transport is indeed reduced 
thanks to the quasi-symmetryii and the turbulent transport 
gains importance. This turbulent transport is thought to be 
driven by microinstabilities like the ion temperature gradient 
mode (ITG) or the trapped-electron mode (TEM). An optimi¬ 
sation to reduce these kinds of microinstabilities in stellarators 
is not trivial, especially since analytical predictions regarding 
the actual nonlinear behaviour of microturbulence in general 
geometry are rather sparse. For instance, only very recently 
Plunk et ah^ published a theory on the saturation of ITG tur¬ 


bulence whereas an analogous theory for TEMs has yet to be 
achieved. Comprehensive numerical simulations of microtur¬ 
bulence in general geometry on the other hand have become 
available, although they are computationally very demanding 
- a typical well-resolved turbulence simulation in flux-tube ge¬ 
ometry and including kinetic electrons needs roughly half a 
million CPUh to reach a saturated state. For this matter, an 
optimisation based on calculating the nonlinear heat flux as a 
figure of merit for every configuration along the path of opti¬ 
misation would be desirable, but it is evident that this proce¬ 
dure is currently not viable. We must therefore And simplified 
expressions to represent the nonlinear heat flux as a figure of 
merit. Ideally, these “proxies” would be based on analytical 
theory of the linear or even nonlinear instabilities. The reduc¬ 
tion of ITG turbulence has been theoretically demonstrated 
via this methodi^^® and in this paper we will tackle the re¬ 
duction of TEM turbulence. Ultimately, one would want to 
combine all methods of optimisation into one grand scheme 
to And the point in configuration space where “the ideal stel¬ 
larator” lives. A code that would be capable of carrying out 
such an ambitious task is STELLOPT'^, which we will also 
use for our TEM optimisation. In the next section we briefly 
explain how STELLOPT works and how a new optimisation is 
implemented. In section |III] we review what we already know 
about TEMs in general geometry using analytical and linear 
numerical findings. Section HVl then shows how we can use 
this knowledge to come up with simple expressions for the 
proxy. In section [V] we present a proof-of-principle equilib¬ 
rium where we optimised starting from HSX towards reduced 
TEM turbulence and we comment on the applicability to ex¬ 
perimentally feasible TEM-optimised equilibria. Section |Vl] 
contains the main conclusions, together with future plans re¬ 
lated to this work. 
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II. OPTIMISING WITH STELLOPT 


III. THE DENSITY-GRADIENT-DRIVEN TEM IN GENERAL 
GEOMETRY 


The STELLOPT code is designed to optimise 3D MHD 
equilibria created by VMEGi^ by minimising the difference 
between certain features of the equilibrium and their tar¬ 
geted value. Each design feature i - this can be the neo¬ 
classical transport, the turbulent transport, ballooning stabil¬ 
ity, the major radius, to name only a few - is associated with 
a target value and the difference between this target 

value and the actual value of the final “optimised” equilibrium 
^equilibrium small as possible. Since usually more 

than only one design feature shall be targeted at once, each de¬ 
sign feature gets assigned a tolerance This acts as a weight 
(1 /ct^) when all design features are eventually combined into 
one function that has to be minimised: 
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( 1 ) 


How ^ tor a given design feature in a given equi- 

librium is determined depends very much on the design fea¬ 
ture itself. Lor the neoclassical transport for example, STEL¬ 
LOPT is coupled to the NEO code^S to calculate the neo¬ 
classical effective ripple ees^. Lor the turbulent transport 
it would be ideal to use the turbulent heat fluxes stemming 
from gyrokinetic simulations. However, since these are CPU¬ 
intensive, simpler proxy functions that can substitute for the 
turbulent heat flux and that are ideally based on analytical the¬ 
ory are sought. The STELLOPT code can be used to opti¬ 
mise any combination of VMEC input parameters to any set 
of target figures of merit (POM), subject to the constraints 
of a given optimisation method. When utilised for stellara- 
tor design the boundary harmonics are treated as the free pa¬ 
rameters, although enclosed toroidal flux, net toroidal cur¬ 
rent and a pressure scaling factor may be included as well. 
As the VMEC boundary representation (R and Z harmonics) 
is non-unique, the initial configuration is converted to either 
Hirshman-Breslau^i or Garabedian22 representation . These 
harmonics are the quantities varied by STELLOPT and con¬ 
verted back to the VMEC representation for evaluation of 
the configuration by VMEC. Once an optimum shape has 
been computed, codes like NESCOIL^^i and C01L0PTi2124 
are utilised to generate a coil set consistent with that equi¬ 
librium. To explore the space of accessible configuration in 
a given device, VMEC may be run in free boundary mode 
and STELLOPT set to treat the vacuum coil currents as free 
parameters. Such a capability allows exploration of a given 
device’s capabilities, as was done for NCSX2^. Addition¬ 
ally, the inclusion of plasma profiles, synthetic MSE diag¬ 
nostics, and magnetic diagnostics^! has allowed the code to 
provide a 3D equilibrium reconstruction capability!^!!!. In 
this paper, a modified Levenberg-Marquardt!! algorithm was 
used to find the minimum of in configuration space. This 
method guarantees that the found optimised equilibrium is at 
least at a local minimum in configuration space. STELLOPT 
is also equipped with stochastic algorithms (e.g. differential 
evolution!!, particle swarm!!) that have not been applied to 
this work. 


To tackle the optimisation of stellarators towards reduced 
TEM turbulence we should first assess what we already know 
about the TEM in general geometry. The TEM!!i!! can be re¬ 
garded as a drift wave that is driven unstable by a resonance 
with the precessional drift of trapped particles. Generally, the 
higher the trapped-particle fraction in a given configuration, 
the more unstable the TEM becomes. It is destabilised by in¬ 
creasing the density gradient and/or the electron temperature 
gradient. In a collisional plasma, trapped particles can become 
detrapped due to collisions, which usually leads to a stabilisa¬ 
tion of the TEM!!. Since this work addresses the suppres¬ 
sion of the worst-case instability, collisions will be neglected 
from here on. Eirst, we revisit some analytical theory regard¬ 
ing the stability properties of TEMs. Because of the temporal 
and spatial scales involved, the gyrokinetic framework is em¬ 
ployed. We then look at different stellarator equilibria with 
very different geometric properties and present linear simu¬ 
lation results that confirm the analytical findings. More de¬ 
tails on the calculations can be found in previously published 
papers!!^—. 


A. Analytical theory 

The stability analysis of TEMs in general geometry via a 
dispersion relation is not as accessible as in tokamaks. How¬ 
ever, it is possible to define a rate of gyrokinetic energy trans¬ 
fer Pe from the fluctuating electric field to the electrons!!^, 
which, at the point of marginal stability where the growth rate 
7 approaches zero, can be written as 


7re^ f dl f _ 

Ue = — — d^U<5(a; -ZUde)wde(wde - W^e)ko</'P/eO- 
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( 2 ) 

Here, Te is the electron temperature, uj is the real frequency 
of the mode and • v^e denotes the bounce-averaged 

precessional drift frequency of the electrons, whose drift ve¬ 
locity is given by v^e- The velocity-dependent diamagnetic 


frequency is given by wj’g = a;*e l + ’7e(^ — |) 
denotes the energy of the particle, rje = ^ d in n 


, where £ 


gives 


dr / dr 

the ratio between the scale lengths of electron temperature 
gradient and density gradient and the diamagnetic frequency 
is defined as (B x kj^) • Vn. In addition, Jq de¬ 

notes the Bessel function, (j> is the electrostatic potential and 
/eo is the Maxwellian distribution function of the electrons. 
Eor the electrons to have a destabilising influence, the en¬ 
ergy transfer rate must be negative, Pe < 0. This means that 
uJdeOJ*e > 0 (if the temperature gradient is small, rje < 2 /3) at 
least for some particles in velocity space, because all the other 
terms are positive definite. Thus, for some of the trapped elec¬ 
trons, the precessional drift must be resonant with the propa¬ 
gation of drift waves. In which direction the trapped electrons 
precess depends on the curvature they sample along their path 
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FIG. 1. The magnetic field strength B (left axis) and the local curva¬ 
ture K (right axis) along a field line of the half-flux flux surface of the 
DIII-D tokamak. The two arrows indicate the paths of two trapped 
electrons with different pitch angle A and thus with different bounce 
points. The deeply trapped particle only samples bad (i.e. negative) 
local curvature along its path, therefore its bounce averaged curva¬ 
ture will also be negative. A barely trapped particle that samples 
both good and bad curvature might have good bounce averaged cur¬ 
vature. More particles are trapped in regions of bad local curvature 
if the two curves of magnetic field and local curvature are in phase. 


in the magnetic field; 


Wde(A) OC 



V'l - ^B{z) 


(3) 


where the integration is taken along the particle path along a 
field line, with the pitch angle like coordinate A = v\/v'^B, 
bounce points Zi and, most importantly, the local radial cur¬ 
vature K, which can have positive and negative values, de¬ 
pending on the direction of the drift. Here, the diamagnetic 
frequency is chosen to be negative, a;*e < 0, which means 
that for a resonance to exist the precessional drift must also be 
negative, ZU^(A) < 0. For a particle to assume such a neg¬ 
ative precessional drift it must sample mainly negative local 
curvature along its path, see Fig. [T] - so-called “bad curva¬ 
ture”. This local “bad curvature” has long been recognised 
as the drive for interchange instabilities^Sili and ITGs, and it 
does play an important role for TEMs, but there its average 
over the bounce motion determines the stability properties. 
Many of the trapped particles will have averaged bad curva¬ 
ture LJd^{X) < 0 if the particles are mainly trapped in regions 
of local bad curvature, i.e. if the magnetic field and the local 
curvature are in phase. Magnetic configurations where this 
is the case should therefore be characterised by destabilising 
electrons, Pe < 0, and should thus be prone to TEM insta¬ 
bilities. Configurations where the magnetic field and the local 
curvature are even partially out of phase, on the other hand, 
should have reduced TEM activity. (Another possibility to 
achieve mainly good average curvature would be to improve 
the local curvature alltogether, for example by having a high 


plasma pressure /3^, but our goal is to also optimise the vac¬ 
uum configurations, so increasing (3 is not an option.) In the 
limit where all particles experience good average curvature, 
uJd^{X) > 0, as is the case in quasi-isodynamic stellarators^^ 
with the maximum-J-property (J is the action integral of the 
bounce motion of trapped particles and constant on flux sur¬ 
faces, the maximum of J being at the plasma centre), it can 
be shown that TEMs and trapped-particle modes are stable in 
large regions of parameter space, i.e. if the electron tempera¬ 
ture gradient is small, pe < 2/3. 


B. Linear simulation results 

Configurations with most of the particles experiencing good 
average curvature can also benefit from enhanced TEM sta¬ 
bility, as can be shown with linear simulations. The simula¬ 
tions are performed with the GENE code^ in the collisionless 
and electrostatic limit. The geometry of the different config¬ 
urations is incorporated into GENE via the GIST geometry 
interface^, and we chose to study three very different stel- 
larator equilibria: the quasi-axisymmetric stellarator design 
NCSX (National Compact Stellarator Experiment, nowadays 
designated QUASAR^-, see Eig. |2] on the left), the quasi- 
helically symmetric stellarator experiment in Madison, Wis¬ 
consin, HSX (Helically Symmetric Experiment, see Eig.|2]in 
the middle), and the stellarator Wendelstein 7-X (W7-X, see 
Fig. |2] on the right), which approaches quasi-isodynamicity. 
For each of the configurations two stellarator-symmetric flux 
tubes were chosen from the flux surface at half toroidal flux, 
s = 0.5, - one where the binormal coordinate a = 0 in the 
midplane, and the second one at a = tt /N where N denotes 
the number of periods. In all configurations, the flux tube 
with a = 0 is centered around the bean-shaped poloidal cross 
section and is therefore referred to as “bean flux tube”. The 
poloidal cross section at the center of the a = ir/N flux tube 
is either triangle-shaped (in W7-X and HSX) or bullet-shaped 
(in NCSX) and therefore called “triangle flux tube” or “bul¬ 
let flux tube”, respectively. The two quasi-symmetric devices 
NCSX and HSX display a strong overlap of the magnetic trap¬ 
ping well and the region of bad local curvature. This is not the 
case though for W7-X, especially at the centre of the flux tube. 
The analytical theory therefore suggests that W7-X should 
have lower TEM growth rates than both NCSX and HSX. On 
the other hand, it should be noted that in NCSX, the region 
of bad local curvature, though it overlaps with the magnetic 
trapping well, is very small compared with the large region of 
good local curvature, especially in the bullet flux tube. NCSX 
might therefore benefit from enhanced TEM stability, too. We 
simulated purely density-gradient-driven TEMs, thus choos¬ 
ing both ion and electron temperature profiles to be flat. Eor 
each value of the normalised density gradient a/L„, where a 
denotes the minor radius and L~^ = — the density gra¬ 
dient scale length, several wave numbers kyPs {ps is the ion 
sound radius) around the expected most unstable mode were 
simulated, and the highest growth rate was then recorded. The 
predicted behaviour of the different configurations is indeed 
born out in the simulations; W7-X and NCSX have the lowest 
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FIG. 2. A comparison of the geometries of the three simulated stellarators NCSX, HSX and W7-X. Displayed are the magnetic field strength 
at the outermost flux surface at the top and the magnetic field strength and local curvature (on the left and right axis, respectively) along the 
two simulated flux tubes per configuration, each of them at half flux s = 0.5. 


linear growth rates, whereas HSX has the highest, see Fig. [3 
The bullet flux tube of NCSX is more stable than the bean flux 
tube, which can be explained by the region of bad local curva¬ 
ture being even smaller in the bullet flux tube than in the bean 
flux tube. The fact that the bean flux tube in HSX has higher 
growth rates than the triangle flux tube can be attributed to 
the fact that there is a magnetic trapping region with bad local 
curvature at zero ballooning angle, which should enhance the 
mode. This is very much in line with the analytical predictions 
and previous linear simulation results^. To summarise this 
section; it was expected from analytical calculations and also 
shown via gyrokinetic simulations that configurations where 
fewer particles have average bad curvature benefit from en¬ 
hanced stability of density-gradient-driven TEMs. Very re¬ 
cent nonlinear results confirm these findings^. More exten¬ 
sive nonlinear data will be published in a later paper. Based on 
this knowledge about the TEM it should be possible to come 
up with a measure, a so-called proxy, that only depends on 
the geometric properties of a configuration and that will give 
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FIG. 3. Linear growth rates of density-gradient-driven TEMs in each 
of the simulated flux tubes in NCSX, HSX and W7-X. At each sim¬ 
ulated density gradient a/L„, where a is the minor radius of the 
device and the density gradient scale length, the growth rate of 
the most unstable mode is displayed. 
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an estimate of the linear (and ultimately also nonlinear) TEM 
activity of the respective configuration. 


IV. PROXY FUNCTION 

The aim of the proxy function is to provide a means of es¬ 
timating the stability of a configuration towards TEMs effi¬ 
ciently, so that the calculation can be performed for many dif¬ 
ferent equilibria during the process of the optimisation. We 
remember the analytical expression for the gyrokinetic en¬ 
ergy transfer rate (Eq. (|2])) and use this as inspiration for our 
proxy function. The central finding from the analytical theory 
discussed above was that it is beneficial for a configuration 
to have as few trapped particles as possible with bad aver¬ 
age curvature. In order to obtain the improved proxy func¬ 
tion Qbounce we thus average the bounce averaged curvature 
7{ oc Wd(A) of a particle with pitch angle A over all trapped 
particles, i.e. over all pitch angles A, equivalent to how the 
average is done in Eq. 0 ; 


Qh 


ounce 






a;d(A)dA, 


(4) 


with 

and where the B^in and B^ax denote the minimum and max¬ 
imum of the magnetic field along a given field line. The mi¬ 
nus sign is introduced to make a configuration with a major¬ 
ity of particles with good average curvature have a minimum 
Qbounce, which seems more intuitive from an optimisation 
point of view. If we plot the maximum growth rate obtained 
from TEM simulations with a pure density gradient for var¬ 
ious configurations and flux tubes versus the corresponding 
proxy value we see that the proxy correlates well with the lin¬ 
ear growth rates, see Eig. |4] Especially if two configurations 
are very different, the proxy correctly predicts which one is 
the more stable. Eor configurations that are very similar, how¬ 
ever, for example different W7-X configurations that mainly 
differ by their mirror ratio (HM being high mirror, EM being 
low mirror, and SC being the standard configuration), a lower 
proxy value does not necessarily mean a lower TEM growth 
rate. This means an optimisation will probably need to make 
large steps in proxy value Q to ensure that the found optimised 
equilibrium indeed has lower levels of TEM activity. 


V. THE PROOF-OF-PRINCIPLE CONFIGURATION 

A first attempt at an optimisation was made with HSX as the 
starting equilibrium. Eor this first proof-of-principle optimi¬ 
sation STELLOPT’s fixed boundary mode was chosen, which 
means the accessible configuration space was very large. This 
was indeed necessary. The constraints of fixed aspect ratio and 
low neoclassical transport prevented STELLOPT from finding 
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FIG. 4. The comparison between the value of the proxy, Qbounce, 
here displayed as Q, and linear growth rate 7 at a density gradient 
a/ Ln = 6 , where the maximum over various binormal wave num¬ 
bers ky was found. Shown are the values for all of the flux tubes dis¬ 
cussed above and additional flux tubes for the low-mirror (LM) and 
standard configuration (SC) of W7-X, where the high mirror (HM) 
is the one that has been used throughout this paper. There is a clear 
correlation between proxy value and linear growth rate between the 
different optimisation families. Note: the high density gradient was 
chosen to ensure the appearance of a strong TEM, especially for the 
fairly stable configurations like W7-X. As can be seen in Fig. [ 3 ] the 
relative stability of different configurations remains the same for the 
higher gradients. 
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Optimised quantity 

Target Count 

Target 

Weight 1/(7 

Neoclassical Transport 

127 

0 

0.001 

Turbulent Transport Qbounce 

25 

0 

1000 

Major Radius Rq 

1 

1.22 

10 


TABLE I. Targets for HSX optimisations. The turbulent and neoclas¬ 
sical values are evaluated at multiple radial locations. 


an equilibrium with lower proxy, i.e. better average curva¬ 
ture. In order to see any change in the proxy the requirement 
of low neoclassical transport had to be relaxed significantly, 
which means that the weight for the neoclassical transport 
was chosen to be very small compared with the weight for our 
proxy (see Table|Il). The resulting TEM-optimised equilibrium 
shown in Eig. |5] on the right has lost the helical symmetry. 
This leads to a significant increase in the neoclassical trans¬ 
port - the neoclassical effective ripple went up by an order of 
magnitude. Moreover, the magnetic field along the field line 
of this preliminary equilibrium is very jagged. This should of 
course be avoided when trying to find a truly optimised con¬ 
figuration. In this case, however, our primary focus is to show 
that the proxy works. To test this, we first performed linear 
GENE simulations. A scan over the binormal wave vector 
kyPs for a purely density-gradient-driven TEM with a density 
gradient a/L„ = 3 and no temperature gradient shows that 
the linear growth rates are indeed reduced for the optimised 
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FIG. 5. Comparison of the initial helically-symmetric HSX and the derived optimised equilibrium produced with STELLOPT. Shown are 
the magnetic field strength on the outermost flux surface at the top and the magnetic field strength and local bad curvature (on the left and 
right axis, respectively) along the bean flux tube at the surface with half flux s = 0.5. The optimised equilibrium is not helically symmetric 
anymore. 


equilibrium, at least for the scales where turbulence is gener¬ 
ated (Fig.|^. This stabilisation also holds for a large range of 
density gradients, as can be seen in Fig. [T] For these simula¬ 
tions, a scan over various wave numbers was performed and 
the highest growth rate for each gradient is displayed. These 
linear results lead to the expectation that a nonlinear simula¬ 
tion of this proxy-optimised configuration would also show 
reduced transport. A test with pure density-gradient-driven 
TEM turbulence at a density gradient of a/Ln = 3 shows that 
the nonlinear electron heat flux went from QIQgb = 1-05 
to QIQgb = 0.62, where the heat fluxes are measured in 
Gyro-Bohm units Qgb = nTf-CsP^s/a?, with the density n, 
the ion sound speed Cs, and the sound Larmor radius ps- This 
means a reduction of about 40% was achieved. However, the 
high neoclassical transport remains a handicap, and, since the 
fixed boundary mode of STELLOPT was chosen to create this 
optimised equilibrium, the magnetic field is not realisable by 
simply adjusting the currents in the existing coils of the HSX 
experiment. Additional STELLOPT runs should therefore be 
used to try to find actually optimised but experimentally real¬ 


isable configurations using the free-boundary mode. 

VI. CONCLUSIONS AND OUTLOOK 

In this paper we have presented a method to optimise stel- 
larators for density-gradient-driven TEM turbulence using the 
optimisation code STELLOPT. The optimisation for TEM tur¬ 
bulence complements ongoing efforts to optimise stellarators 
not only for neoclassical transport but also for turbulent trans¬ 
port. We used analytical theory and linear flux-tube simula¬ 
tions performed with the GENE code to guide us in devising a 
proxy function that can stand in for the expected turbulent heat 
flux. Both analytical theory and the linear simulations sug¬ 
gested that configurations with a lower fraction of particles 
with bounce-averaged bad curvature should be less unstable 
to TEMs. The bounce-averaged curvature averaged over all 
trapped particles was thus chosen for the proxy. The compar¬ 
ison between the proxy and linear growth rates for TEMs in 
various configurations revealed that the proxy is well suited 
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TEM with L„=3.0, L,. =0.0 




FIG. 6. A comparison of the linear TEM growth rates 7 and the real frequencies cj in the bean flux tubes of the initial HSX equilibrium and 
the optimised equilibrium for different wave numbers kyps at a fixed density gradient a/Ln = 3. The optimised equilibrium has lower growth 
rates at the turbulence relevant scales. 


TEM a/L„ only in the initial and optimised HSX configuration 



FIG. 7. A comparison of the linear TEM growth rates in the bean flux 
tubes of the initial HSX equilibrium and the optimised equilibrium 
at different density gradients ajLn- The optimised equilibrium has 
reduced growth rates for all density gradients. 


to predict the relative stability of a configuration. Assuming 
that the linear growth rates are correlated with the turbulent 
transport levels, the proxy should thus be able to guide the 
optimiser towards configurations with lower TEM turbulence 
levels. A first proof-of-principle configuration where this was 
indeed achieved was presented. There, the linear growth rates 
were reduced compared with the starting equilibrium of HSX, 
as was the turbulent heat. This configuration was, however, 
not realisable with HSX’s given coil set. The presentation 
of an experimentally feasible TEM-optimised configuration is 
deferred to a future publication. 



f(a/Cj) 


FIG. 8. A comparison of the nonlinear heat flux in the bean flux 
tubes of the initial HSX equilibrium and the optimised equilibrium 
at a density gradient a/L„ = 3. The grid for these simulations 
was chosen as follows: in the radial direction nkx = 192, in the 
binormal direction nky — 48 with the lowest wave number being 
kymin = 0.01, along the field line nzO = 64 for the initial config¬ 
uration or nzO — 256 for the optimised configuration, respectively, 
for the parallel velocity nvO = 40 and for the magnetic moment 
nwO = 20. 


One possible improvement of the current proxy would be to 
include further weighting of the deeply trapped particles by 
taking into accound the mode structure of the linear modes 
via the electrostatic potential |(/)p, as it is included also in the 
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equation for the energy transfer rate, Eq. (|2l). Another pos¬ 
sibility that has proven fruitful in the reduction of ITG tur¬ 
bulence is to include the distance of the flux surfaces in the 
optimisation, trying to find configurations where this distance 
is particularly large, which would result in a smaller effective 
gradient and thus possibly large regions of stability in parame¬ 
ter space. Combining the TEM optimisation with the ITG op¬ 
timisation where the local curvature is minimised might also 
lead to better results for TEM turbulence. The simultaneous 
reduction of both ITGs (and interchange instabilities in gen¬ 
eral) and TEMs might already be happening with our proxy, 
when the reduction of the bounce averaged bad curvature is 
achieved by making the local curvature better. However, there 
might be configurations where this is not feasible, but shifting 
the local bad curvature away from the magnetic wells is. 

In addition to improving the proxy function and thus hopefully 
reducing the turbulence at high gradients other problems in 
turbulence optimisation could and should be addressed. One 
of these problems is to increase the critical gradient for the 
onset of turbulence, which is particularly important if the tur¬ 
bulence is very “stiff”, i.e. if the heat flux increases dramat¬ 
ically once the critical gradient is exceeded. Ideally, in ad¬ 
dition to reducing the turbulent heat flux, one would achieve 
an increase in the particle flux to flush out impurities from 
the plasma, but a deeper understanding of turbulence is re¬ 
quired before this challenge can be tackled. In general it re¬ 
mains to be seen to what extent the simultaneous optimisation 
of several different aspects (neoclassical and turbulent trans¬ 
port, fast-particle confinement, divertor etc.) can be success¬ 
ful. Until then the TEM optimisation presented here will serve 
as a useful tool to learn more about the influence of geometry 
on TEM stability and will certainly guide future efforts. 
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